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ABSTRACT 

Context. Globular clusters show star-to-star abundance variations for light elements that are not yet well understood. The preferred 
explanation involves a self-enrichment scenario, within which two subsequent generations of stars co-exist in globular clusters. 
Observations of chemical abundances in the main sequence and sub-giant branch stars allow us to investigate the signature of this 
chemically processed material without the complicating effects caused by stellar evolution and internal mixing. 
Aims. Our main goal is to investigate the carbon-nitrogen anti-correlation with low-resolution spectroscopy of 20-50 stars fainter than 
the first dredge-up in seven Galactic globular clusters (NGC 288, NGC 1851, NGC 5927, NGC 6352, NGC 6388, and Pal 12) with 
different properties. We complemented our observations with 47 Tuc archival data, with four additional clusters from the literature 
(M 15, M 22, M 55, NGC 362), and with additional literature data on NGC 288. 

Methods. In this first paper, we measured the strengh of the CN and CH band indices, which correlate with the N and C abundances, 
and we investigated the anti-correlation and bimodality of these indices. We compared icn, the ratio of stars belonging to the CN- 
strong and weak groups, with 15 different cluster parameters. 

Results. We clearly see bimodal anti-correlation of the CH and CN band stregths in the metal-rich clusters (Pal 12, 47 Tuc, NGC 6352, 
NGC 5927). Only M 15 among the metal-poor clusters shows a clearly bimodal anti-correlation. We found weak correlations (sligthly 
above 1 cr) of icn with the cluster orbital parameters, present-day total mass, cluster concentration, and age. 

Conclusions. Our findings support the self-enrichment scenario, and suggest that the occurrence of more than two major generations 
of stars in a GGC should be rare. Small additional generations (< 10-20% of the total) would be difficult to detect with our samples. 
The first generation, which corresponds to the CN-weak stars, usually contains more stars than the second one (<rcw>=0.82±0.29), 
as opposed to results based on the Na-0 anti-correlations. 

Key words, stars: abundances - stars: main sequence -globular clusters: individual (M 15, M 22, M 55, 47 Tuc, NGC 288, NGC 362, 
NGC 1851, NGC 5927, NGC 6352, NGC 6388, NGC 6752, and Pal 12) 



1. Introduction 

Recent work revealed that several Galactic globular clusters 
(hereafter GGC) had a much more complex star-formation his- 
tory than previously thought. This has implications for their use 
as simple stellar population (SSP) templates to study more com- 
plex and distant stellar systems and galaxies, and therefore re- 
quires careful study. 

Based on deep and accurate photometric studies, GGC were 
found to contain different stellar populations, in the form of 
multiple evolutionary sequences in their color-magnitude dia- 
grams (CMD). The first clus ter found t o hos t multiple stel- 
lar populations was to Cen (Lee et al.', '1999!; 'Pancino et al.L 
2000 ; Bedin et al., 2004; Sol lima et al., 2005; Vil lanova et aL, 
20071). Lately NGC 2808 dPiotto et alfTioOTr NGC 185 1 
dMilone et all 120081; IVentura et al.L 120091; iHan etaLl l2009h . 



NG C 6388 (iMoretti et al.L l2009l). M 22 dPiottol l2009h . and 47 
Tuc dAndersijn et all 20091) joined the rapidly growing group of 
complex stellar population GGC. 

But even before the multiple photometric sequences were 
detected, anti-correlations between the strength of the CN and 
CH bands around 3880, 4200, an d 4300 A were found, start- 
ing with the pioneering work of lOsbornI ( 1197 ll) . in all prop- 
erly observed clusters. These anti-correlations are sometimes 
bimodal in nature (see Sect. 14.11 for references), with two sep- 
arate groups of CN-strong and CN-weak stars. The first de- 
tections w ere of c ourse l imited to the red giant branch (RGB) 
stars (see ISmithL I1987L iKraftL 11994 for early reviews), be- 
cause of the instrumental limitations at that time, but later 
they were also found among sub giant br anch (SGB) st a rs and 
even among main sequence (MS) stars ( Cannon et all 1998 : 



* Based on FORS observations collected at the European Southern 
Observatory, Chile, within the observing programs 68.D-0510 and 
69.D-0056. Also based on data obtained from the ESO Archive, within 
the observing program 67.D-0153. 



CoSenL iHarbeck. Smith & GrebelL 120031: ICarretteetal 



20051) . Anti-correlations were found among other light ele- 
ments such as Na, O, Al, Mg (see the revi ews bv iKraft. 1994 ; 
Gratton. Sneden & CarrettaL |2004 and F dCunha et all I2OO3I; 
Smith et all l2005h . 
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Table 1. Observing logs (see Sect. |2] and |5]for details). 



Cluster 


[Fe/H] 
(uex) 


(m-M)v 
(mag) 


E(B-V) 
(mag) 


run 


^ Stars 




(sec) 


S/N(3800A) 


NGC 104 (47 Tuc) 


-0.68 


13.37 


0.04 


2001-07 (Archive) 


64 


3 


11700 


20.5±7.8 


NGC 288 


-1.39 


14.83 


0.03 


2001-10 


72 


5 


19800 


23.6±4.8 


NGC 1851 


-1.22 


15.47 


0.02 


2001-10 


50 


3 


18000 


28.7+3.6 


NGC 5927 


-0.55 


15.81 


0.45 


2002-07 


45 


4 


32400 


9.9+1.9 


NGC 6352 


-0.54 


14.44 


0.21 


2002-07 


24 


1 


9900 


19.0+3.0 


NGC 6388 


-0.44 


16.14 


0.37 


2002-07 


22 


2 


19800 


8.9+2.5 


NGC 6752 


-1.50 


13.13 


0.04 


2001-10 


56 


4 


7200 


28.0+5.7 


Pal 12 


-0.80 


16.47 


0.02 


2001-10 


23 


2 


9900 


14.1+3.9 



No sample of Milky Way field stars observed so far showed 
any sign of anti-correlation s (but they do s how s i gns of dredged 
up CNO processing, see 



; (but tney go s now s i gns or aredged 
Hanson et all ll998UGratton et al.L 



120001 iMishenina etall l2006l) . Also, no chemi cal anomalies 
were f ound in open cluster s (see references in Ide Silva et al.L 
120091: iPancino et al.l l2009l) or in the field stars belonging 
to dwarf galaxies in the Local Group. Only recently, Na- 
O anti-correlations were detected in extragalactic GC (glob- 
ular clusters) of the Fornax dwarf galaxy and of the Large 
Magellanic Clotids dLetarte et all l2006t iJohnson et all \2QQd 
iMucciarelh et al.L I2009I) . This suggests that environment plays 
a fundamental role in the occurrence of these anomalies. 

The presence of CH-CN anti-correlations was soon inter- 
preted as a signature of CN(0) cycle processing, which tends to 
deplete C and to enhance N, and of mixing that can bring CN(0) 
processed material to the stellar surface, during the first dredge- 
up that happens at the base of the RGB. To reach this conclu- 
sion, it is as sumed that CN tra ces the N abundace and CH the C 
abundance dSmith et al.Lll996h . This is well supported by abun- 
dance analysis work based on spectral synthesis, which revealed 
[C/Fe1 and [N/Fe] spr eads of up to mor e than 1 dex (see e.g., 
iBrilev et all r2004a.b; Cohen et all l2005h . However, this intrin- 
sic scenario, often referred to as the stellar evolution scenario, 
was soon put into trouble by two observational facts: (i) CH and 
CN anti-coiTelations in MS stars cannot be produced by the stars 
themselves, because they have low mass (M^O.8 M©) in GGC 
and do not burn hydrogen through the CNO cycle; ( it) the pres- 
ence of Na, Al, Mg, and O (anti-)correlations points towards the 
NeNa and MgAl chains that take place at much higher temper- 
atures ([Penisenkov & Denisenkova, 1990; Langeretal., 19931 
iPrantzos et al.L 2007h . which are not reached in the MS or RGB 
phases of GGC stars. The addition of inter nal pollution mecha- 
nisms such as rotation-induced mixing (e.g. lSweigart & MengelL 

—I I J ( ; 

1979ljCharbonnel, 1995) or the so-called canonical extra-mixing 



(e.g. lDeni ssenkov & VandenBerg, 2003i) did not solve the prob- 
lem because they still could not explain the Na-O and Al-Mg 
anti-correlations. Of course, mixing does o c cur among red giants 
above the base of the RGB (ISmith et al.L I1996L iGratton et all 
I2000L ISmith & MartellL l20()l . and this is the reason why it is 
better to focus on stars on the MS and SGB, which still have not 
undergone any mixing episode, to have a clean perspective on 
chemical anomalies, independent from stellar evolution effects. 

A natural alternative was then some sort of extrinsic sce- 
nario, where the chemical anomalies are produced outside the 
G GC stars. An exampl e is the primordial scenario proposed 
by ICannon et al ] (Il998h . which foresaw the presence of two 
clouds with dififerent chemical composition, merging in the for- 
mation phase of GGC. Another example is the stellar pollu- 



tion scenario dP'Antona et al l ll983LlThoul etal.Ll2002L) . where 
only the stellar surfaces are polluted by CN(0) processed winds. 
These and other extrinsic s cenari os ar e extensively discuss ed by 
iHarbeck. Smith & Grebell (l2003h . and lKavser et"al] (I2008L) . 

The most promising of the extrinsic scenarios remains 
the so-called self- e nrichm ent scenario, first proposed by 
ICottrell&Da Costal (1198 Ih . which assumes the presence of two 
(or more) subsequent stellar generations, the first polluting the 
intra-cluster medium from which the second one forms. A very 
important constraint to this scenario is that most GGC show a re- 
markably homogeneous composition as far as heavier elements 
(e.g., Ca, Fe) are concerned. In the past five years, the chase has 
thus moved from finding the most promising scenario to seeking 
for the first generation stars responsible for polluting the intra- 
cluster medium. Three proposals have been made so far: ( i) in- 
termedi ate mass (M>3-5 Mp) asymptotic giant branch ( AGB) 
stars (e.g. Ventura et al., '200 iL iDenissenkov & HerwigL 12003 



Ventura & D'Antona, 2009), which naturally explain chemical 
anomalies, but also provide clear links to the helium problem, 
proposed to explain the multiple populations found in some 
clusters (see above) and their horizontal branch (HB) mor- 
phologies ( see the series of pape rs by Carretta et al., the last 
ones being Carretta et all 2009allbl); (ii) fast rotating massive 



stars dMaeder & Mevnetri2006L iPrantzos & CharbonnelL I2006L 
iDecressin et al.L 2007 ). which suffer from meridional circula- 
tion, able to bring CN(0) processed material to the surface, and 
to release it through slow winds, through their accretion disc and, 
finally, through explosion; this scenario helps also in solving 
the helium and multiple population problems; a nd ( Hi) the very 
recently proposed massive interactive binaries (Ide Mink et al.L 
2009), which are known to lose mass and could indeed explain 
the observed anomalies. Moreover, the slow winds of all these 
types of polluters are prone to diff'erent degrees of stripping de- 
pending on the cluster and its environment. This can explain 
why different clusters show different behaviours of their mul- 
tiple populations and chemical anomalies, as described in detail 
in Sect. |6] 

The aim of this paper is to look for CH and CN band strength 
variations and bimodalities in a sample of unevolved stars in 12 
GGC with different properties. Those clusters that show clear 
variations and bimodalities will be further studied in a following 
paper dedicated to abundance determinations, to investigate the 
extent of the [C/Fe] and [N/Fe] variations. Main sequence stars 
have not undergone any dredge-up phenomenon yet, so they are 
free from any evolutionary effect that, superposed to the CH and 
CN extrinsic variations, could confuse the picture. We will also 
correlate the CH and CN strength variations with GGC param- 
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Fig. 1. Colour-magnitude diagrams for the newly observed GGC in our sample. Grey dots show our pre-imaging V, V-I photometry, 
caUbrated with l i teratu re data (see text). For 47 Tuc, which was not part of our original observations, we used the photometry by 
[Rosenberg et alj (120001) . A selection on the photometric parameters from DAOPHOT II is applied to show only the best measured 
stars. Black dots mark spectroscopic targets. Clusters are ordered as a function of their metallicity from the most metal-poor (top 
left) to the most metal-rich (bottom right). 



eters, in order to provide additional observational constraints to 
the self-enrichment scenario. 

The plan of the paper is the following: in Sect. |2] we discuss 
the observations and data reductions; in Sect. [3] we define the 
indices passbands and their measurement; in Sect. |4] we present 
our results; in Sect.|5]we discuss the results on a cluster by clus- 
ter basis; in Sect. |6] we correlate the CH and CN strengths with 
cluster parameters, and in Sect.|2]we summarize our findings and 
draw our conclusions. 



2. Observations and data reduction 

We selected a sample of eight GGC with different properties to 
test the correlations of their CH and CN strengths with clus- 
ter parameters and environmental conditions; 47 Tuc (from the 
ESQ archive), NGC 1851, NGC 288, NGC 5927, NGC 6352, 
NGC 6388, NGC 6752, and Palomar 12. We also tried to select 
stars covering a large area in each cluster, to be free from biasses 
due to possible ra dial trends. Measure ments for four more GGC 
were added from'Kavs er et al.1 (l2008h : M 15, M 22, M 55, and 
NGC 362, plus measurements for additional stars in NGC 288 
(see Sect. I3.3l l. A more detailed discussion of each cluster's 
properties and results can be found in Sect.|5] 



2.1. Observations 

Observations were carried out mostly in service mode in two 
different runs and two pre-imaging runs at the ESO VLT/UT2 
telescope in Cerro Paranal, Chile, with the FORS2 multi-object 
spectrograph dAppenzeller et al in MXU mode. The ob- 

serving logs can be found in Table [1] The first run was dedi- 
cated to the three metal-poor clusters (NGC 6752, NGC 288, 
and NGC 1851) and to Pal 12, and was carried out between 
September and October 200 1 , with relatively good conditions for 
the pre-imaging (sky mostly clear and seeing ^1") and less fa- 
vorable conditions for the spectroscopy observations (sky mostly 
veiled and seeing ranging from 0.5" to 2"). The second run was 
dedicated to the metal-rich clusters (NGC 6352, NGC 5927, and 
NGC 6388) and was carried out between April and July 2002. 
The conditions for both pre-imaging and spectroscopy obser- 
vations where less fortunate than the in the first run, with the 
sky often veiled or cloudy, and seeing between 1.5"-2". This 
explains in part why the colour-magnitude diagrams (CMD) of 
these clusters are of a much worse quality (Fig.[I]) and why the 
S/N ratios of their spectra tend to be quite low as well (Table [T]i. 

Although we used the same grism, GRJS_600B, the old 
2kx2k SITE detector (24 yum/pix) was replaced by a mosaic 
composed by two 2kx4k MIT CCDs (15 //m/pix) between the 
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two runs. This implies that the covered wavelength ranges were 
not exactly the same in the two runs; also, the red sensitivity 
of the CCD was improved, but this had no effect on our mea- 
surements. The typical resolution of our spectra in both runs is 

2.2. Pre-imaging reductions 

We obtained a set of images for one or more fields within each 
cluster, generally consisting of a short exposure (a few seconds) 
and a longer one (a few tens of seconds) for each of the Johnson 
V and the Cousins I filters. All images were bias-subtracted and 
flat-fielded with the ESO FORS pipeline. Point spread function 
(PSF) fitting photometry was carried out with the DAOPHOT II 
ALLSTAR and ALLFRAME packages (IStetsonl [19871 [19941) 
using a constant model PSF across the field, which experience 
showed to yield the best results. 

The photometric calibration was done using stars in common 
with photom etry from the lite r ature: we used the V, I catalogue 
pubHshed by Bellazzini et al.' ('2001') for NGC 1851; the B, V 
catalogue by .Dalessandro et al., C2008b) supplemented by I mag- 
nitudes (Dalessandro, priv ate communication) for NGC 6388; 
and the V, I catalogue by [Rosenberg et al.l (l2000l) for 47 Tuc 
(NGC 104). For all rernaining clusters we used the HST pho- 
tometries bv lPiotto et alJ (l2002h . The resulting calibrated colour- 
magnitude diagrams and programme stars are shown in Fig.[T] 



2.3. Spectroscopy reductions 

The spectroscopic target stars were selected from the pre- 
imaging data as the most isolated stars located around the turn- 
otf (TO) and SGB regions of the CMD (Fig. [T]i- In most cases, 
more than one exposure was taken with a single mask, to reach 
a higher S/N ratio. With each mask, we were able to observe 
20-30 stars with 1" wide, 10" long slitlets. The average S/N for 
each cluster (after the frame pre-reductions) in the CN 3880 A 
region, together with the total number of good stars analysed in 
each cluster (see Sect. 12.51 ), is indicated in the last two columns 
of Tabled] 

For the data pre-reduction, we used IRAlQ for overscan cor- 
rection and bias-subtraction. We removed cosmic ra y hits with 
the IRAF Laplacian edge-detection routine dvan Dokk un, 
The frames were subsequently flat-fielded and reduced to one- 
dimension spectra with the IRAF task apall. The sky background 
was subtracted using the information on both sides of the stellar 
spectrum. Finally, for those fields that were observed more than 
once, we co-added all spectra together on a star by star basis. We 
checked that the shift between the different spectra (< 1 km s ' ) 
was negligible compared to our wavelength calibration uncer- 
tainty (^40 km s"') before coadding. 

The shape of the final spectra is affected by the different 
response of the instrument at different wavelengths. There are 
two ways to remove this effect: flux calibration and contin- 
uum normalization. The first needs observations of standard 
stars, which were not available in our case. The alternative is 
to fit a polynomial to the pseudo-continuum and use it to nor- 
malize the spectrum. However, the CN 3880 A band region 
contains many atomic lines and molecular bands, and the re- 
sponse of the instrument is quite low in this ultraviolet re- 
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Fig. 2. Influence of telescope and/or instrument flexures on the 
derived radial velocities. Each point represents one star, with its 
radial velocity and the X and Y coordinates of the corresponding 
slitlet center. The solid Unes represent a bilineal fit to the data 
(see text). 



gion. Therefore, no attempt was made to normalize the spec- 
tra, as m any other authors skip this step for the same rea- 
son (e.g. ICohen et al.L 12002; H arbeck. Smifli & GrebelL [20031: 
[Kavser et al.L [20081 among others). 



' IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 



2.4. Radial velocity and shifts 

Finally, radial velocities were calculated by cross-correlation 
(with fxcor in IRAF), choosing the highest S/N star on each 
MXU plate as a template. The templates Vr were computed us- 
ing the laboratory positions of the strongest lines (e.g. H^, H^, Ca 
H and K) in our wavelength region with the IRAF task rvidlines. 

The three dominant sources of uncertainty for the final V,- 
of our spectra are: (i) uncertainty in the wavelength calibration; 
(ii) small shifts in the star centering within each slitlet; and (iii) 
variations across the field of view owing to instrument flexures. 
The uncertainties in the wavelength calibration are generally re- 
moved using the positions of strong emission sky lines, which 
could also be used to evaluate th e influence of the instrument 
flexures (e.g. iGallart et al.L [200 ll) . Unfortunately, in the wave- 
length region observed here, there are no strong sky lines that can 
be used. The uncertainty due to small misalignment of stars in 
the slitlets could be estimated with the through-slit images (see 
ICarrera et aU[2007[) . and we found it negligible (<10 km s ') in 
comparison with the Vr errors (~20^0km s"', see below). 

Significant variations of Vr across the field of view, most 
probably due to flexures, are clearly observed as shown in Fig. [2] 
The estimated radial velocity varies systematically with the po- 
sition of each slitlet along the Y CCD axis (and on the X axis as 
well), by as much as ~80 km s"', i.e., larger than the the velocity 
uncentainty (~20-40 km s"') and of the scatter expected within 
each cluster (of the order of a few km s '). The example in Fig. [2] 
shows a typical difference of about ^60 km s To correct for 
these variations, we computed a bilineal fit for each observed 
plate, in the form Vr-a + bX + cY (see Fig.|2]i. Stars near the cen- 
tre of each MXU mask, which have the smallest shifts, were used 
to roughly r eport all the m ask velocities to the GGC systemic ve- 
locity (from [Harrisl[l996[) . The final velocities corrected for the 
above variations are reported in Table |2| 



2.5. Quality control and sample selection 

Given the huge amount of spectra (see Table [T]) and the some- 
times low S/N achieved during observations, we appUed the fol- 
lowing selection criteria: 
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Table 2. Index measurements and radial velocities for the sample stars. A portion of the table is shown for guidance about its 
content, the complete table is available in electronic format through the ADS service. 



Cluster 


Star 


RA 


Dec 


CN 


6CN 


eivcN 


CH 


(5CH 


errcH 


Vr 








(deg) 


(deg) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(kms-') 


(kms-') 


NGC 104 


0101 


6.58283 


-71.88794 


-0.262 


-0.266 


0.029 


0.837 


0.872 


0.040 


-52 


33 


NGC 104 


0102 


6.58437 


-71.88708 


-0.313 


-0.236 


0.035 


1.024 


0.953 


0.058 


-23 


40 


NGC 104 


0104 


6.59971 


-71.88342 


-0.290 


-0.147 


0.059 


1.020 


1.016 


0.085 


-40 


52 


NGC 104 


0105 


6.63296 


-71.88597 


-0.404 


-0.128 


0.061 


1.077 


1.024 


0.106 


-36 


49 


NGC 104 


0106 


6.65038 


-71.88544 


0.074 


-0.184 


0.060 


0.978 


0.987 


0.069 


-44 


45 


NGC 104 


0107 


6.58404 


-71.87014 


-0.454 


-0.310 


0.018 


0.902 


0.877 


0.042 


-12 


28 


NGC 104 


0109 


6.64762 


-71.87739 


-0.393 


-0.247 


0.033 


1.000 


0.938 


0.064 


-19 


31 



- spectra with S/N<8 (per pixel) in the CN 3880 A band region 
were rejected; 

- spectra with significant defects (spikes, holes) in the mea- 
surement windows were rejected; 

- stars with discrepant radial velocity by more than 2.5 cr from 
the cluster average were rejected; 

- stars with discrepant Ca(H-i-K) index measurement were re- 
jected (see Sect. [3]for details); 

- stars with discrepant Hg index measurements were rejected 
(see Sect. [3]for details); 

The criterion that had the highest impact on star rejections 
was by far the first one. The final sample contains 356 stars (out 
of an initial set of ~600) belonging to eight clusters, as indicated 
in the third column of Table [T| and also reported in Table |2l 



3. Index definition and measurement 

The strength of a spectral feature can be evaluated with a spec- 
tral index, which is defined through one window centred on the 
atomic line or molecular band that we will evaluate, and one or 
more windows around it, used to define the continuum level. In 
the literature we can find many definitions of the spectral indices 
used to measure the strength of the CN and CH molecular bands. 
In most cases, these definitions are optimized for red giant stars, 
observed at a specific spectral resolution. In our case, because 
we observe subgiants and dwarfs, and we aim to compare our 
results wit h others in the literature, we chose to adopt the indices 
defined bv lHarbeck. Smith & Grebel (l2003h : 



S3839{CN) = -2.5 log 
54142(CA^) = -2.5 log 
Ci/4300 = -2.5 log 



^3861-3884 
^^3894-3910 

^4120-4216 

0.5/^4055-4080 + 0.5/^4240-4280 

/^4285-4315 

0.5/^4240-4280 + 0.5/^4390-4460 




4100 4aoo 
Wavelength (A) 





3900 4000 

Wavelength (A) 



4B40 4B60 4BB0 
Wavelength (A) 



Fig. 3. Two spectra of CN-weak, CH-strong stars with different 
S/N ratio are shown as examples of our best and worst data: 
the upper spectrum in all panels has S/N^25 (on the S3839(CN) 
band) and belongs to NGC 1851; the lower spectrum has S/H^ 1 
and belongs to NGC 6388. Both are shifted on an arbitrary flux 
scale for sake of clarity. The top panel illustrates the chosen win- 
dows used for the CN and CH indices (magenta hatched regions) 
along with their respective continuum windows (grey hatched re- 
gions); while the bottom panels show the windows adopted for 
the H and K Calcium index (bottom left) and the H^ index (bot- 
tom right). 



where /'3861-3884, for example, is the summed spectral flux 
in ADU counts from 3861 to 3884 A. The related uncer- 
tainties have been obtained with the expression derived by 
IVolImann & Eversb erg (200^, assuming pure photon (Poisson) 
noise statistics in the flux measurement^ Fig- E] (top panel) 

- We could neglect here the readout noise, always between 3 and 5 
ADU, compared to our 300-3000 ADU gathered in the CN band region; 
the error in the sky subtraction, with tens of pixels on each side of the 
spectrum for sky estimation, and an order of 100 ADU sky level for the 



shows the adopted indices windows; examples of our best and 
worst S/N spectra are overplotted. 

We also measured two additional indices, centred around the 
calcium H and K lines and the H^ line (see also Fig. [3] bottom 
panels): 



worst cases; and the spectrum tracing errors, because we always had 
a peak signal - in the central pixel - of at least 70 ADU, and a total 
300-3000 ADU level, even in the CN region that had the lowest S/N. 
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HK = 1 



^39 10-4020 

^^4020-4130 



^^4847. 



875-4876.625 



0.5/^482 



7.875-4847.875 



+ 0.5/^4876.625-4891.625 



The H and K line strengths depend mostly on temperature and 
on the calcium abundance, while the line strength depends on 
temperature and of course on the hydrogen abundance. We used 
these two indices to reject a few remaining 3 cr outliers (at most 
2-3 stars per cluster). Since we used the indices only in a rela- 
tive sense, we did not need to compare our measurements with 
the literature, therefore we defined our own narrow and conser- 
vative windows (see Fig. |3] lower panels) in order to minimize 
the spread due to disturbing features containing other elements 
and to maximize our ability to pinpoint outliers. 

3.1. CN-band choice 

In our wavelength region, we measured two different CN band 
indices, S3839 for the CN band aro und 3880A and S4142 for 
the m uch weaker one around 4200 A jHarbeck. Smith & Grebell 
l2003h . Fig. El shows a comparison between the two. The cor- 
relation appears weak for the metal-poor clusters (NGC 6752, 
NGC 288, NGC 1851), in spite of the higher S/N of the spectra, 
because these double-metal molecular bands are weaker A clear 
correlation appears for the most metal-rich clusters (47 Tuc, 
NGC 6352, NGC 5927, and NGC 6388), although it is not as 
striking for NGC 6388, given the low S/N and paucity of stars. 
The slope of the relation between S4142 and S3839 appears to 
be roughly the same for all clusters (dotted line in Fig.|4]i, while 
the zeropoint varies slightly from cluster to cluster 

As expected, Fig.|4]clearly shows that S4142 is less sensitive 
to the CN variations, because generally its spread is just slightly 
larger than the median errorbar shown in the lower-right corner 
of each panel. This is especially true for the metal-poor clusters. 
On the other hand, the spread of S3839 is always larger than its 
median errobar, even for the metal-poor clusters. For this rea- 
son, we will be relying on S3839 measurements for CN, and we 
will set aside the less reliab le S4142 o nes (as done also by, e.g., 
ICannon et ailil998;,Harbeck. Smith & Grebel .20031) . 



3.2. Dependency on temperature and gravity 

At a fixed overall abundance, it is well known that the CN- 
band and CH-band are stronger in stars with lower tempera- 
ture and gravity. To eliminate this dependency, different au- 
thors use diffe rent proxies for temperature and gravity, such 
as c olours (e.g. Harbeck. Smith & Grebel, 2003) or magnitudes 
(e.g. lNorris et alT 19811) or so me combination of the Balmer in- 
dices (e.g. lKavser et al.Ll2006h . The indices are then corrected by 
fitting the lower envelope of the distribution in the chosen plane, 
and the new indices are generally indicated as <5S3839(CN) and 
5CH4300. Most of the cited authors find generally no significant 
trend for S3839(CN), but they do always correct the CH4300 
index, usually with a second-order polynomial. Given the diver- 
sity of shapes and slopes, we prefer to use median ridge line^to 



^ A rough estimate of the uncertainty in the placement of the me- 
dian ridge line can be obtained by using the first interquartile range of 
the rectified indices, divided by the square root of the total number of 
points. We find a typical uncertainty ranging approximately from 0.01 
to 0.02 in the CN index, and from 0.005 to 0.01 in the CH index, de- 
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Fig. 4. Comparison between the CN S3839 and the S4142 in- 
dices. Each panel shows measurements for the indicated cluster, 
and the median errobar is also shown in the lower right corner. 
Clusters are sorted by metallicity, from the most metal-poor (top- 
left) to the most metal-rich (bottom-right). Dotted lines are lin- 
ear fits with a fixed slope (that of the whole sample) and rougly 
adjusted zeropoints. 



correct for the curvature induced by both temperature and grav- 
ity effects. We show an example of the computed median ridge 
lines, with NGC 1851 in Fig. 13 

Generally, we confirm that the curvature of S3839(CN) as 
a function of V magnitude is much smaller than the curvature 
of CH4300. We define our rectified 5S3839(CN) and 5CH4300 
indices as the difference between the original S3839(CN) and 
CH4300 indices and their respective median ridge lines. The rec- 
tified indices will be used throughout the rest of this paper 

3.3. Additional data from tine literature and arcliives 

To increase our observed sample, we re-reduced the 47 Tuc data 
we found in the ESO archive with our procedure, which is very 
similar to the one adopted by the original tea m that made the 
observations dHarbeck. Smith & GrebelL l2003h . This way, a di- 
rect comparison of the results was possible, showing excellent 
agreement and therefore providing a useful validation of our 
adopted procedure. We also added index measurements for the 
four clusters in Kayser et al. (2008) that had MS and SGB stars 
faint enough to have pristine surface composition, not altered by 



pending on the cluster. These uncertainties are negligible for all appli- 
cations in the present work. 
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Fig. 7. Anti-correlation plots for the CH and CN band strengths. Each panel shows the measurements for stars (grey dots) in each 
cluster CH strong and weak stars are separted by the horizontal dotted line, and their centroids are marked as large white dots, while 
CN strong and weak stars are separated by the vertical dotted line and their centroids are marked as large black dots. The S/N ratio 
in the CN band region is indicated in the lower left corners. 



the first dredge-up mixing episode, which happens at the base of 
the RGB. We therefore selected from their measurements (their 
Table A.l) those stars in NGC 362 with V<18.5, in M 55 with 
V<17.5, in M 15 with V< 1 8.5, and in M 22 with V<17 (see 
Fig. 2 bv lKavser et all l2008l) . Their definition of indices is ex- 
actly the same as ours (Sect. |3]l. We applied a more strict se- 
lection in radial velocity, at 2.5 cr (as in Sect. 12.5b . We then re- 
moved the temperature and gravity dependency from their mea- 
surements (their Table A.l), using the same procedure described 
in Sect. [32] Kayseretal. (2008) do not give the S/N of their 
spectra, but judging from their errobars (see their Table A.l), 
they should be of the order of S/N^15-20 for M 15, M 22, and 
M 55, and S/N:^ 10-15 for NGC 362 anf NGC 288. In the fol- 
lowing sections, we analyse and discuss these data together with 
ours. 



4. Results 

4.1. CH and CN bimodalities 

The presence of CH and CN bimodalities, i.e., two well sepa- 
rated groups of stars with a different strength of these bands, 
was apparent almost as soon as CH and CN anomalies were 



discovered in red giants in the 1970s and 1980s. Some clus- 
ters showed clear bimodalities in CN (or in [N/Fe]) and 
weaker or missing bimodalities in CH (or in [C/Fe]), while 
some other clusters instead showed a continuous spread in CN 
strength. Various studies i dentified CN bimodalities among re d 



2004b), 



giants of M 3 and M 13 dSuntzeffl [1981'; 'Bril ev etal. 

NGC 6752 (iNorris et al.L Il981h. u j Cen (.Cohen & Bell 11986 ). 
NGC 6934 (ISmith& Belli 1 19861) NGC 6171 (ISmithL ll988h■ 



M71 (ISmith & Penny. '1989^. 'Lee. '2005'). M 2 (ISmith & Mated 
FT990I) . M 5 (Ramirez & Cohen, 2002}), NGC 28 8 an d 
NGC 362 jKavser et al.i . 2008), NGC612j_(Marino et al.Ll2008l) . 
NGC 6356, and NGC 6528 (Martell & Smith..2009.) . among oth- 
ers. Bimodalities among MS and sub- giant bra nch stars have 
also been found in M 71 ( Cohenl [19991) . 47 Tuc ( Cannon et al 



1 1998: Harbeck. Smith & Grebek I2003L _ 

ICarretta et al.Ll2005h . and NGC 6752 dCarretta et aUl2005h 



Brilev et all l2004a : 



The interest related to the presence of a bimodal distribu- 
tion - rather than a continuous spread - is clear in the light of 
the latest theories of self-enrichment for GGC (see also Sect.[T|i. 
First generation stars should have a "normal" surface compo- 
sition, rich in C and poor in N, while the second generation 
ones, formed from CN(0) c ycle processing pollute d gas, should 
have higher N and lower C (ISmith & Norri sL ll982h . If we find a 
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Fig. 5. Example of the removal of temperature and gravity de- 
pendencies from the CH and CN indices, using median ridge 
lines on NGC 1851. Large black dots represent the target stars, 
while small grey dots represent the calibrated pre-maging pho- 
tometry of NGC 1851. The median ridge lines of CH and CN as 
functions of V magnitude are marked as dotted lines. 




S/N{CN)=15-eO 



S/N(CN)=30 




_ S/N{CN)=10-25 



S/N(CN)=:30 




(!S3839(CN) 



(!S3B39(CN) 



Fig. 8. Zoomed version of those clusters that show a small spread 
on the scale of Fig.|2] 



clear bimodality in the data, this supports the idea of two stel- 
lar generation^ If a bimodal distribution is found measuring 
the CH and CN band strengths, we expect that the underlying 
[C/Fe] and [N/Fe] distributions will turn out to be bimodal as 
well, since each cluster has the same overall metallicity and all 
stars in each cluster have roughly the same atmospheric parame- 
ters. This of course requires confirmation via spectrum synthesis 
calculations. These calculations will be the subject of a future 
paper Indeed, such a bimodality in C a nd N enhancements has 
been found e.g., bylBriley et alJ (l2004bt) in RGB stars in M 13, 
or bv lBrilevet aD(l2004ah in MS stars in 47 Tuc. That such bi- 
modalities have rarely been found among red giants when the 
Na-O anti-correlations are considered deserves further studies. 
We point out two possible reasons in advance: (i) the Na and es- 
pecially the O abundances in giants are difficult t o measure, and 
indee d we often see only upper limits for [O/Fe] (ICarretta et al.L 
l2009aiibi) or (ii) C, N, O, and Na are not directly comparable, 
because they are not produced exactly at the same temperature: 
C and N are altered within the CN bi-cycle, O is depleted in the 
complete CNO cycle and Na is produced in the NeNa cycle, each 
dominating at progressively higher temperatures; in this second 
hypothesis much care should be taken with the comparison of 
Na-O and C-N anti-correlations with each other. Finally we note 
that if such a C-N bimodality is present, it should be easier to 
see in MS stars than in RGB stars, where it could be attenuated 
by internal mixing (Sect.lTJ. 

We studied the presence of bimodal distributions for the CN 
and CH rectified band strengths by means of histograms (Fig.|6]l 
and - more importantly - by plotting our measurements in the 
CH-CN plane (Fig. [T] see also next Sect.). We consider a dis- 
tribution to be "bimodal" when the centroids of the CN-strong 
(CH-weak) and CN-weak (CH-strong) stars are clearly sepa- 



The opposite is not necessarily true, in the sense that self- 
enrichment could also produce a continuous spread in the CH and CN 
band strengths, depending on the details of the chemical enrichment 
history, such as duration and intensity of star-formation episodes. 



ratec0 in the CH-CN plane and - as a supporting evidence - 
when the histograms also show signs of bimodality. Detailed 
studies of the statistical significance of such bimodalities are de- 
ferred to a following paper, dealing with the [C/Fe] and [N/Fe] 
distributions. First, we note that all the most metal-rich clusters 
(Pal 12, 47 Tuc, NGC 6352, and NGC 5927), with the only ex- 
ception of NGC 6388 (which has low S/N and few stars), show a 
clear bimodality not only in ^S3839(CN), but also in 5CH4300. 
Among the metal-poor clusters, the situation is less clear There 
is a bimodality only in M 15, even if it is the most metal-poor of 
the sample and has a relatively low S/N ratio. For the remaining 
clusters we only see vague hints of possible bimodalities (see 
Sect.|5]for more details) in Fig.|6] but not in Fig.|7]or|8] We will 
discuss individual clusters in Sect. |5] 

4.2. CN and CH anticorrelations 

While the S/N of our metal-poor clusters spectra (NGC 362, 
NGC 6752, NGC 288, NGC 1851, M 22, M 55) is compara- 
ble to that of other a uthors (i.e.. lHarbeck. Smith & Grebeill2003l: 
iKavser et al.Ll2008h . the double molecule bands of CN are really 
weak at those low metallicities, for MS stars. On the contrary, 
even if the S/N of our remaining metal-rich clusters is lower due 
to bad weather conditions, CH and CN bands are much stronger. 
For the MS stars extracted from the lKavser et al.l (l2008h sample, 
the S/N is relatively low because they focussed their analysis on 
brighter RGB and SGB stars. 

Because our S/N is sometimes at the limit of detection in 
the CN band region, at least for the most metal-poor clusters, 
we used the 5CH4300 index, which has the highest S/N, to 
separate target stars in a CH-strong (5CH4300>0) and a CH- 
weak (5CH4300<0) group for each cluster A weighted aver- 
age was used to compute the centroids of the CH-strong and 



' This generally implies a significant centroids separation of at least 
1 cr, but this separation of centroids does not grant a bimodality per se, 
as for NGC 6388. 
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CH-weak groups in the CN-CH plane. The resulting centroids 
with their 1 cr errors are reported in Fig. [T] as large white dots 
along with measurements for each star. For completeness we 
also divided the stars in CN-strong (5S3839(CN)>0) and CN- 
weak ((5S3839(CN)<0) groups, and their centroids are marked as 
large black dots with their errobars. A zoom-in for those clusters 
with small spreads on the scale of Fig.|2]is presented in Fig. [8] 

We see clear anti-correlations (Fig.|7ll only for those clusters 
that showed a clear bimodality in both indices, i.e., M 15, Pal 12, 
47 Tuc, NGC 6352, and NGC 5927. We would be tempted to 
say that where our data are good enough to show an anticorrela- 
tion, this anticorrelation is always clearly bimodal. The case of 
NGC 288 seems to be an exception, and it is further discussed 
in Sect.|5] In most of the other cases, anti-correlations and/or bi- 
modalities were found by other authors, as discussed in detail in 
Sect.|5j so we suspect that in general the present data do not have 
high enough S/N to be conclusive on these clusters. 



5. Cluster by cluster discussion 

5.1. NCG 104 (47Tucanae) 

The cluster NGC 104 or 47 Tuc is one of the most massive and 
meta l-rich GGC (Table [3]l . It is considered to be part of the thick 
disc jDinescu et al with a low eccentricity and relatively 

high inclination orbit. Currently it is located ^ 3 kpc below 
the Galactic plane in project ion on the sky close to the Small 
Magellanic Cloud. iPritzl e^ al. ( 200 5|) place i t in th e group of 
thick-disc clusters, while Mackev & Gilmord (l2004l) classify it 
as bulge/disc cluster. 

It is one of the most studied GGC thanks to its vicinity and 
its relatively high [Fe/H]. While the abundance determinations 
on red giants are too man y to mention here (the most recent be- 
ing lCarretta et aUl2009cl) . the abundance analysis of turn-off and 
early subgiant s tars re sulted i n rFe/H1^-() . 67 ±0. 01(±0.04) dex 
(ICarretta et al.L |2004 . while iJames et al.l (|2004 derive -0.69 
±0.06 dex for SGB stars and -0.68+0.01 dex from three turnoff 
stars. It was one of the first GGC in which CN-CH vari- 
ations were detected among MS stars (ICannon et al 1 [19981 
[firilev et al., 1994). The star-to-star scatter has been clearly de- 
tected for light elements; while the Na-O anticor relation is not 
very s trong, it has nevertheless been detected by ICarretta et al.l 
(l200l . 

The CN-CH molecular abundances of turn-off and sub-giant 
stars in this cluster were also topic of a detail e d stud y by 
iHarbeck. Smith & Grebell (l2003h and iBrilev et al.l (l2004ah . re- 
sulting in a clear bimodality in band strengths and also in [N/Fe], 
with abundance differences of up to 1 dex. We use here their 
archive data to test the homogeneity of reductions. A clearly bi- 
modal anti-co rrelation is detected here (Fig.|7 1, wh ich confirms 
the resu lts bv lHarbeck. Smith & Grebell (120031) and lBrilev et al.l 
(l2004al) and which, together with the double SGB and wide MS 
found by Anderson et al. (2009), demonstrates the presence of 
two distinct populations in this cluster 



5.2. NGC 288 

The source NGC 288 is a low concentration, low central den- 
sity cluster located close to the South Galactic Pole, in a very 
low interstellar extinction region. Its retrograde motion and rela- 
tively low velocity com ponents place it amon g the inner old halo 
(OH) globular clusters (iDinescu et al. I IT997I) . Tidal shocks due 
to bulge and disc for this cluster are estimated to be strong com- 



pared to the usually more predom inant internal relaxation and 
evaporation (IDinescu et al.Ll 19991) . 

Its horizontal branch is composed almost entirely of stars 
bluer than the instability strip. It is often studied together with 
NGC 362, which has the same metallicity, but very different HB 
morphology: they represent one of the classical second param- 
eter pairs. £olte (1989J was one of the first to propose helium, 
CNO, or age differe nces as possible explanatio ns for the second 
parameter problem. IShetrone & Keanel dlOOOh made a compar- 
ative spectroscopic study of red giants in the second-parameter 
GC pair NGC 288 and NGC 362. The [Fe/H] they derived for 
NGC 288 is -1.39 + 0.01 dex. These authors also remark on 
the Na-O and Al- O anti-correla t ions, w hich were found in both 
clusters. Fig. 7bv lCarretta et alj ( l2009bh shows even a hint of bi- 
modality in the Na-O anti-correlation. A clear bimodality in the 
C H and CN index stre ngth has been detected am ong red giants 
bv lKavser et~al] (l2008h and lSmith & Langland-Shulai(t2009i) . 



Our data are the ones by iKavser et al.l (l2008l) . of which we 
selected only unevolved stars, and we do not see any clear anti- 
correlation. There is no reason why unevolved stars in a GGC 
should have a uniform composition while the evolved ones show 
a bimodality, because no known mixing mechanism would be 
able to create a bimodality from a homogeneous population. 
Therefore we must conclude that our data do not have a suffi- 
cient S/N to reveal a bimodality in these metal-poor MS stars. 
Spectral synthesis calculation will be able to quantify the sensi- 
tivity of our spectra to C and N abundances, while more higher 
quality data are most probably needed anyway to settle the ques- 
tion. We note here that as we will see in Sect. |6l NGC 288 falls 
perfectly in line with other clusters in all correlations with clus- 
ter parameters. 



5.3. NGC 362 

As mentioned before, NGC 288 and NGC 362 are one of the 
second-paramete rs pairs , with similar metallicity and different 
HB morphology. iBoltd (Il989l) proposed CNO, among others, 
as one of the second p arameters. NGC 362 has a metallic- 
ity of -1.33+0.01 dex (IShetrone & Keanel l2000l) . it is more 
concentrated than NGC 288 (Table [3]) and slightly closer to 
the Galactic centre. The anomalies in the CH and CN band 
strengths among its giants were studied more t han in NGC 288 
(McClure & Norrisl 119741: iFrogel et all 119831; ISmith & NorrR 
liSJ ISmithi Il98l . iKavser et al.l (l2008l) extensivelv studied 
both NGC 288 and NGC 362, finding a clear anti-correlation 
among RGB and SGB stars in both clusters, and a clear bimodal- 
ity among red giants in both clusters. This suggests that anti- 
correlations or chemical anomalies cannot be the predominant 
source of the different HB morphology in these two clusters. 



As for NGC 288, bimodalities and an anti-correlation were 
found in the CH and CN band strengths of RGB stars, while our 
data for the MS stars reveal none. We suspect that here also this 
is simply owing to the low S/N of the spectra (lower than for 
NGC 288), and it must be worsened by the smaller sample. A 
study with higher quality spectra would be even more interesting 
than for NGC 288, because NGC 362 does not behave like the 
other GGC in the correlations with cluster parameters (Sect. ^ 
and it is not clear whether NGC 362 is a real outlier or if the 
present data are simply inconclusive. 
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(dex) 
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(Gyr) 
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(M0L3') 


NGC 104 (47 Tuc) 


-0.75" 


-0.99 


3.756' 


13.06 


0.09 


11.5" 


2.03 


42.86 




-9.42 


6.05'i 


1.331 


NGC 288 


-1.26" 


0.98 


4.22 l-i 


10.62 


0.09'" 


2.9" 


0.96 


12.94 


0.093 


-6.74 


4.851 


2.151 


NGC 362 


-1.33'' 


-0.87'' 


4.079' 


10.37 


0.01 


6.4" 


1.94 


16.11 


0.023 


-8.41 


5.53'i 


0.90' 


NGC 1851 


-1.22''' 


-0.32'' 


4.097' 


9.98 


0.05 


10.4" 


2.32 


11.70 


0.015 


-8.33 


5.491 


1.611 


NGC 5927 


-0.67'' 


-1.00 


3.724' 


12.67 


0.04 


4.3° 


1.60 


16.68 




-7.80 


5.32' 




NGC 6352 


-0.55= 


-1.00 




12.67 


0.07 


5.4° 


1.10 


10.51 




-6.48 


4.57' 




NGC 6388 


-0.41" 


-0.69'' 


4.255' 


12.03 


0.01 


18.9" 


1.70 


6.21 




-9.42 


6.021 


1.891 


NGC 6656 (M 22) 


-1.76' 


0.91'' 




12.67 


0.14 


9.0" 


1.31 


28.97 




-8.50 


5.561 


2.071 


NGC 6752 


-1.56" 


1.00 


4.47 P 


11.78 


0.04 


4.5" 


2.50 


55.34 




-7.73 


5.16' 




NGC 6809 (M 55) 


-1.97" 


0.87 


4.153) 


12.29 


0.02 


4.9" 


0.76 


16.28 




-7.55 


4.991 


3.231 


NGC 7078 (M 15) 


-2.34" 


0.67'' 


4.477' 


12.93 


0.05 


12.0" 


2.50 


21.50 


0.027 


-9.17 


5.84' 




Pal 12 


-0.80*5 


-1.00 




8.83 




1.3° 


1.94 


17.42 


0.193 


-4.48 


3.751 





NOTES: All parameters are derived from the 2003 revision of the 'Harris' ("1996) GGC catalogue except: " 'Carretta et al.' ('2009c'); ^ 
[Shetrone & Keane (2000); = Yona et al. (2009); Kraft & Ivans (2003); = Feltzing et al. (2009); ' Marino et al. (2009); » Cohen (2004); '' Cateli3 
120091) ; '|R ecio-Blanco et alj fe006l) ; ' B ragaglia (2009 , priva te c ommunication); ' relativ e a ges by Marin-Fr anch et alj j2009l) conver t ed to ab- 
solute ages multiplying by 12.8 Gy r; "' [ Frenk & Falll (Il982j) : " iPrvor & Mevlaiil ( Il993l) ; ° Ibnedin et al. (2002h ; p iBaumgardt et all ( l2009h ; 1 
iMcLaughlin & van der Marell ilOOS) : 'iMandushev et al. ( 199]]). 

Meaning of columns. (1) cluster name (alternate name in parenthesis); (2) mean iron abundance ratio; (3) HB morphology, where HBR=(B- 
R)/(B-hV-i-R), being B the number of stars bluer than the instability strip, R redder, and V the number of variables in the strip; (4) the maximum 
temperature of the HB; (5) the age in Gyr; (6) the isophotal ellipticity 6=l-(b/a); (7) the central radial velocity dispersion; (8) the concentration 
(2.5 means core collapsed cluster); (9) the tidal radius; (10) the ratio between the half-light and the Jacobi radius; (11) the integrated V magnitude; 
(12) the logarithm of the total mass in solar units; (12) the mass-to-light radius in solar units. 



5.4. NGC 1851 

The cluster NGC 1851 has a very eccentric and highly in- 
clined o rbit, with a large excursion into the outer parts of the 
Galaxy dPinescu et al.Lll997l Il999l) . The iron abundance mea- 
surement s from high-resol ution spectra rep ort rFe/H1=-l.22 + 
0.03 dex (lYong et al.Ll2009h and -1 .27+0.09 (lYong & Grundahli 
l2008h . The bimodal HB in this cluster already hi nted towards 
multip le stellar populations, but only recently iMilone et al.l 
( I2OO8I) detected two parallel S GB sequ ences, thanks to exquisite 
ACS@HST photometry, and iHan et al. (2009) found that the 
use of U, U-I photometry reveals a striking bimodality in the 
RGB. A poss i ble di fferent radial distribution of the two SGB 
(IZoccali et al.Ll2009l) would support some pre dictio ns of the self- 
enhchment scenario (iD'Ercole et al.L l2008h . but iMilone et al.l 
(2009) do not confirm this finding. 

A possible explanation for the double SGB is a 
CNO enriched stellar co m ponen t (ICassisi et al 1 I2OO8I: 
ISalaris. Cassisi & Pietrinferril l2008h . possibly - but not 
necessarily - combined with some age difference of the order 
of 1 Gyr at most. This is expected to produce a strong effect 
on CN and CH molecules. Indeed, among the R GB stars some 
extrem ely CN strong stars have been detected (iHesser et al.L 
Il982h . Another study reporting a large C NO abundance spread 
in giant stars is that bv lYong et al.l (l2009l) . The HB morphology 
of NGC 1851 seems to rule out a large He difference be tween 
the two populations (ISalaris. Cassisi & PietrinferniLl2008h . 

Unfortunately, as for NGC 288 and NGC 362, at this metal- 
licity our S/N ratio is probably not sufficient to reveal anti- 
correlations and bimodalities. We only see some bimodality in 
the 5CH4300 index (the one with the highest S/N), but not in 
5S3839(CN), and only in Fig.|6] but not in Fig.|7]or||] Moreover, 
our targets do not fall in the ACS field, so we cannot attempt any 
experiment on their position on the two different SGB found by 
iMilone et al. at this stage. More studies including higher 

S/N spectra and a careful target selection could help in clarifying 



Table 4. Clusters sample positional and orbital parameters. 



Name 


Class"-^ 


Rgc 








P 






(kpc) 


(lO^km^s--) 


(kpc) 


(deg) 


(10«yr) 


NGC 104 


BD (TK) 


7.4= 


-872= 


3.1= 


29= 


190= 


NGC 288 


OH (H) 


12.0= 


-787= 


5.8= 


44= 


224= 


NGC 362 


YH (H) 


9.4= 


-856= 


2.1= 


21= 


208= 


NGC 1851 


OH(-) 


16.7= 


-340= 


7.6= 


22= 


584= 


NGC 5927 


BD (~) 


4.5= 


-1020' 


0.7' 


9' 


147' 


NGC 6352 


BD (TN) 


3.3= 










NGC 6388 


BD (-) 


3.2= 










NGC 6656 


OH (TK) 


4.9= 


-871= 


1.9= 


18= 


190= 


NGC 6752 


OH (TK) 


5.2= 


-977= 


1.6= 


18= 


156= 


NGC 6809 


OH (H) 


3.9= 


-1038= 


3.7= 


56= 


122= 


NGC 7078 


YH(H) 


10.4= 


-752= 


4.9= 


36= 


242= 


Pal 12 


SG(H) 


16.2'' 


650'' 


20.1'' 


58'' 


730'' 



NOTES: Literature ref e rences are: " GGC classification by 
iMackev & van den BergH ( |2005[) . where BD=bulge/disk, OH=old 
halo, YH=youn g halo, SG=Sagi ttarius; '' GGC classification (in 
paranthesis) by iPritzl et al.l ( 120051) , whe re TK=thick disk. TN=thin 
disk, H=halo ; = 20 03 revision of the iHarrisI ( Il996l) catalogue; '' 
iDinescu etal] (I2OOOI); = result s with model JSH95 by iDinescu etaTI 
( I1999I) : ' ICasetti-Dinescu et al.l ( l2007h 

Meaning of columns. (1) cluster name; (2) classification according 
to Galactic sub-component; (3) Galactocentric radius; (4) total orbit 
energy; (5) maximum height above the Galactic plane; (6) orbit 
inclination; (7) orbital period. 



the nature of these substructures. As for NGC 288, NGC 1851 
behaves like all other clusters of similar metallicity when trends 
with GGC parameters are considered (Sect.|6]l. 



12 



Pancino et al: CN and CH band strengths of MS stars in 12 GGC 



5.5. NGC 5927 

The source NGC 5927 belongs to the metal-rich globular clus- 
ters. It suffers from relatively high reddening and thus was not 
well st udied spectros c opical ly with high-resolution data: accord- 
ing to IZinn & Wesll (11984 . it is one of the most metal-rich 
disc clusters ([Fe/H]=-0.3 dex): ICohenl (Il983l) found this clus- 
ter to b e -1-0.59 dex more metal rich than 47 Tuc: lKraft & IvansI 
derive from low-resolution spectra [Fe/H]//=-0.67 dex; 
[Francois (1991) derived [Fe/H] -1.08 dex based on the high- 
resolution study of one star only. Given the large spread in litera- 
ture determinations, we adopt the lKraft& IvansI (l2003h estimate, 
keeping in mind that NGC 5927 could well be more metal-rich 
than that. To the best of our knowledge, the present study is the 
first one dedicated to light element anti-correlations in this clus- 
ter 

We find a quite clear and bimodal anti-correlation of the CH 
and CN band strength, in spite of the low S/N ratio of our spectra 
and of the relatively high reddening in the cluster field. Clearly, 
at this metallicity the CN and CH bands are so strong that it 
is quite easy to reveal these variations. Finding a bimodal anti- 
correlation in such a metal-rich cluster is interesting, since for 
exa mple no open cluster revealed these variations up to now 
/seelde Silva et al.L 120091: fPancino et al.Ll2009HMartell & SmithL 



I2009r and references therein). Obviously, high-resolution spec- 
troscopy of stars in this neglected cluster would be of enormous 
interest. 

5.6. NGC 6352 

NGC 6352 is a sparsely popul ated metal-rich cluster with disc 

ace it in th e group of thin 



kinematics. iPritzl et al 



i copul ated 
20051) pli 



disc clusters, while Mackev & Gilmor^ (l2004l) classify it as a 
bulge/disc GG C. Based on high-resolution spectra of eight stars, 
ICohenldTgsl found this cluster to be -1-0.38 dex more metal rich 
than 47 Tuc. iCarretta & Gratton ( 199 7j) re-analysed equivalent 
width measurements of iGrattonI Il987h for three stars finding 
an aver age [Fe/H]= -0.64 + 0.06, and recentlv [Carretta et alj 
(l2009d) re-adjuste d it to [Fe/H]^ -0.62 + 0.05. A recent high- 
resolution study bv lFeltzing et al.l (l2009l) . based on nine horizon- 
tal branch stars, found [Fe/H]= -0.55 + 0.03 and suggested that 
on the HB there is a correlated trend in Al and Na abundances, 
similarly to what was found for stars on the RGB in other glob- 
ular clusters. 

No study dedicated to li ght elements was ca rried out to 
our knowledge, so the hint bv iFeltzing et al.l (l2009l) is our only 
comparison. We do find a clear and bimodal CH and CN anti- 
correlation, which is even more striking if one considers the 
low number of stars observed and the relatively low S/N around 
3900 A. As for NGC 5927, the high metallicity helps in reveal- 
ing CH and CN band strength variations, and the study of these 
metal-rich GGC would be interesting for the same reasons. An 
HST photometric study - seeking for multiple sequences - cou- 
pled with high-resolution spectroscopy for deriving the chemical 
patterns of NGC 6352 could be one of the next steps in the study 
of chemical abundance anomalies in GGC. 



5.7. NGC 6388 

The cluster NGC 6388 is one of 10 most massive in the Milky 
Way, it is a bulge cluster, located about 3 kpc from the Galactic 
centre. It is very centrally concentrated and tightly bound, and 
has amon g the highest predicted escape velo city at the clus- 
ter center dMcLaughlin & van der MarelLl2005h . Based on high- 



resolution FLAMES -UVES spectra. lCarretta et al.l(l2007bh mea- 
sured [Fe/H]= -0.44 + 0.01(±0.03). These authors also detected 
the presence of Na-O and Mg-Al anti-correlations, and of a Na- 
Al correlation among RGB stars. A cl ear bimodality in CH and 
CN w as also detected on the RGB by ISmith & Langland-Shulal 
(l2009l) . 

NGC 6388 is s ometim es referred to as being unusual (e.g. 
iDalessandro et al.l, l2008bh : in contrast to expectations for its 
high me tallicity, the cluster harbours an extended blue horizontal 
branch (iRich. Minniti & Lieberl Il993h . and RR Lyr stars with 
periods much longer than expected for their metallicity. In addi- 
tion, the horizontal branch presents a slope, so that in the V-band 
its blue tail lies a bout 0.5 mag brighter than the red HB clump 
(iRaimondo et al.L 12002) . These features are not reproducible by 
stellar evolutionary models with an SSP of standard GGC abun- 
dance ratios, but could be explained with the presence of more 
than one stellar pop ulation and some s elf-enrichment with CNO 
processed material (lYoon et al.L 12008"). The evidence for multi- 
ple stellar pop ulations came from t he most recent deep near-lR 
photometry of iMoretti et al.l ( l2009l) . who detected two distinct 
sub-giant branches in this cluster. 

Unfortunately, the strong differential reddening and the in- 
sufficient S/N of our data prevent us from reaching any conclu- 
sion about the CH and CN correlation at the MS level. But given 
its high metallicity and total mass, this is one of the most inter- 
esting clusters for further investigations. In spite of the poor data 
for this cluster, NGC 6388 follows all the trends of other GGC 
with clusters parameters (Sect.|6]l. 



5.8. M 22 (NGC 6656) 

Besides o) Centauri, M 22 is one of the first GGC that 
were suspected to host mult i ple ste llar populations and chem- 
ical anomalies dHesser et al.L Il977l) . Unfortunately, str ong dif- 
feren tial reddening has complicated its study (Richter et al.L 
1 19991) . therefore spectroscopic studies initially gave conflict- 
ing results: while a spread in the CH a nd CN abundances 
of R GB stars appeared unquestionable dNorris & FreemanL 
119831) . some studies reported on abundance variations of 
0.3+0.5 dex in Ca and/or Fe, often correlated with the CH 
and CN vari a tions (TPetersonL ll980L iPilachowski et all I1982L 
iLehnert et all |1991L iBrown & WallersteinL 1 19921) . while other 
studies found no significant variation i n the heavy e l ement 
content jM anduca & BellL I1978L ICohenL ll98lL iGrattonL I1982L 
iLaird et al.h J991i) . This was probably because the reported vari- 
ations were of the o rder of th e quoted uncertaintie s (see also 
iMonaco et all |2004 . Recently, iMarino et all (l2009l) found that 
two groups of stars exist in M 22, where [Fe/H] appears corre- 
lated with the ^-process elements and wit h [Ca/Fe], and the usua l 
Na, O, Al anti-correlations were found. iDa Costa et all (l2009l) 
confirm that the Ca distribution shows a large spread with a pos- 
sible bimodality. 



Our data, belonging to the unevolved stars from the sam- 
ple by iKayser et al. {200^ do not have sufficient S/N to reveal 
CH and CN anti- correlations o r bimo dalities. No clear bimodal- 
ity was found bv lKavser et al.l (l2008l) among RGB stars, either. 
If really different metallicities co-exist in the cluster, this could 
further complicate the detection of anti-correlations and of their 
bimodalities, if they are present. Most probably spectral synthe- 
sis and the determination of [C/Fe] and [N/Fe] would help in 
clarifying the picture. 
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5.9. NGC 6752 

The cluster NGC 6752 is one of the best studied regarding its 
chemical properties, given its relative proximity to us (Table [1]). 
It has an intermediate metallicity a nd a blue HB. Its orb it is re- 
markably similar to that of 47 Tuc (iDinescu et all 12000'). and it 
is therefore associated with the thick disc by some authors. Many 
high-resolution spectroscopic studies of this cluster provide the 
following iron abundance measurements: [Fe/H]= -1.48±0.01± 
0.06 dex by (Gr atton et al ., 200^ based on FLAMES-UVES 
spectra of 7 giants near the RGB bump, [Fe/H] = - 1 .42 based 
on U VES spectra of MS and sub-giant stars (iGratton et al.L 
I2OOI . [Fe/H]= -1.49 ± 0.07 for 9 SGB stars, and [Fe/H] = 
-1.48±0.07 for9 tu rnoff stars based on FLAMES-UVES spectra 
(iJames etal.U2004 . [Fe/H]= -1.61 based on 38 RGB stars ob- 
served at very high resolution dYong et al.Ll2005l) . and [Fe/H]/ = 
-1.56, [Fe/H]// = -1.48 based on FLAMES G IRAFFE spectra 
of 137 RGB stars dCarretta et al.Ll2007ai 120093) . 

The O-Na anticorrelation is well estab lished in NGC 6752 
from observat i ons o f both MS and SGB ( iGratton et all 12001 ; 
Carretta et all l2005h. a s well as of RGB (iNorris & Da Costa , 



1995HYo nget al.'.l 2003h . Data from iGrundahl et al.l (l2002h and 
Yong et al . (2003) show Na-O and Al-Mg anticorrelations in 



stars that are both b righter and fainter than the RGB bump. 
iPasquini et al.l (l2008h detected very strong N abundance varia- 
tions i n the MS Turn-off stars. Furhtermore. lYong & Grundahll 
(l2008h detected for the first time correlations of light ele- 
ments with Si and hea vier elements among red giants. Recently, 
IVillanova et al.l ( 2009h have obtained a direct He measurement 
from the 5875 A line in HB stars and found a homogeneous 
abundance of Y=0. 245+0. 012, thus putting in doubt He as the 
origin of multiple populations, at least in this GGC. 

Unfortunately, our data have certainly too low a S/N to 
clearly detect anti-correlations. The histograms in Fig. |6] are 
suggestive of a possible bimodality, but we know akeady 
that the ant i -corre lations are bimodal in this cluster because 
INorris et"aL l(fT98l foun d a clear bimodality among RGB stars 
and lCarretta et all (l2005h among MS stars. 



5.10. M 55 (NGC 6809) 

The source M 55 is a moderately sized cluster (Table [3]), with 
an unusually low central concentration, which makes it rela- 
tively easy to collect stellar samples even within the cluster core. 
The cluster lie s in front of the southern tail of the Sagittarius 
dwarf galaxy (iLane et al I I2OO9I) . There are few spectroscopic 
meas urements of metalUcity for M 55 stars, yiel ding fFe/H1=- 
1.65 dCaldwell & DickensLll988h . [Fe/H] =-1.95 dMinniti et al.L 
Il993h . and recently [Fe/H]=-1.93 +0.02(±0.07) dex fro m 14 
stars observed with UVES@VLT (ICarretta et al.L l2009d) . who 
also found a Na-O a nti-correlation and some ind ication of pos- 
sible Al variations. arin-Franch et al. I (I2OO9I) assign M 55 
to the old group of GGC, based on their relative ages scale. 
The blue straggler stars (BSS) of M 55 suggest an unusually 
high primordial popu lation of binaries compared to other GGC 



dLanzoni et al.L 2007h. although this contrasts with the conclu- 
sions of SoUima et al.l d2007h Another interesting result con- 



ce rns the abundance of helium , meas ured with the R-method 
bv I Vargas Alvarez & San dquisti d2007h . which appears to be the 
highest recorded for a massive GGC: Y=0.274+0.016. 

Our results on this cluster are somewhat unclear, with a pos- 
sible secondary peak in the CH band histograms, but not in the 
CN ones (Fig. (6]). This is most probably owing to the insuffi- 
cient S/N ratio for MS stars in this metal-poor cluster. However, 



0.4 



o.a - 





-1.5 

[Fe/H] (dex) 



-0.5 



Fig. 9. Relation between the CN-strong and CN-weak centroids 
geometrical distance in the CN-CH plane and [Fe/H] from 
Table[3] The white dot represents M 15, which is excluded from 
the linear fit (dotted line) of the remaining clusters (black dots). 



iKavser et al.l d2008l) found no clear sign of anti-correlation, even 
consid ering the hig h er S/N red giants spectra, and an early pa- 
per by iBrilev et"al] dl993h suggests that the CN spread might 
be unusually low in this cluster. So the question of CH and CN 
anti-correlations and of their bimodalitiy in M 55 remains open. 



5.11. M 15 (NGC 7078) 

The most metal-poor cluster in our sample is M 15, with 
[Fe/H]=-2.34 (Carretta et al., 2009c,), it has a blue HB and it 
is a core-collapsed GG C, which is suspected to harbour a black 



IS a core-collapsed UU C, which is suspected to harbour a blacij 
hole in its centre (see iHo et al.L I2003L iMcNamara et alL 12003 



van den Bosch et a U I2006L IChakrabartvl l2006t 1 Kiselev et al 



2008r and references therein). It is at present quite far away from 
the Galactic center (Table|4li. 

Anti-correlatio ns were detected am ong all usual light ele- 
ments in M 15 (see lSneden et all 1 19971 and references therein), 
with hints that there could b e a sp read also in the heavy ele- 
ments (see also lSneden et al.L I2OO O: 'Otsuki et al.'.'2006'). On the 
low -resolution spectroscopy front. .Cohen et aL (i2005i) obtained 
a beautiful C-N anti-correlation among stars on the lower RGB, 
with impressive abundance variations of 1-2 dex. However, 
neither their [C/Fe] nor [N/F e] distributions show any evi- 
dence of bimodality. Also. lLed d2000i) found an anti-correlation 
am ong bright gia n ts. The Na-O anti-correlation was also found 
bv ICarretta etal] d2009bl) . although their measurements were 
mostly upper limits. 

We selected our tar gets among the unevolved stars measured 
bv lKavser et al.l d2008l) . who found some bimodality among both 
evolved and unevolved stars. With our median ridge line correc- 
tion of the temperature and gravity effect (Sect. 13.21 ) we were 
able to reveal a clearly bimodal CH and CN anti-correlation. 
This is really striking if one considers the very low metallicity 
and the S/N ratio of these spectra, lower than that of NGC 1851 
and comparable to that of NGC 288. Indeed, the index strengths 
measured in this paper are - as expected - a function of global 
metallicity, represented by [Fe/H] in Fig. A roughly linear 



* ACN is the geometrical centroids distance of the CN-strong and 
CN-weak groups in the CH-CN plane. The geometrical distance is 



14 



Pancino et al: CN and CH band strengths of MS stars in 12 GGC 



trend appears when the distance between the centroids of the 
CH-weak and CN-strong stars is plotted against [Fe/H], with 
some sign of "saturation" at the metal-poor end, where the in- 
dices become weak. The behaviour of M 15 clearly stands out 
of this simple relation, and a prediction of the present work is 
that its [C/Fe] vs. [N/Fe] anti-correlation should be significantly 
more extended than those of other clusters in the sample. 

5.12. Palomar12 

A very interesting cluster to include in our sample is Pal 12, 
because it is one of the GGC traditionally associated with 
the Sagittarius dwarf spheroidal galaxy (Sgr), which is at 
present disrupting under the strain of the Milky Way tidal field 
hbata et al.', 1995; Maiewski et al.','2003!). Its lower [a/Fe] ratio 
(iBrown et al.,.1997: Coheru,2004i,) and younger age than typical 



Table 5. Ratio of CN-strong and CN-weak stars and centroids 
distance. 



GGC dRosenberg et al. , 19981: iMarm-Franchet al its lo- 

cation close in Galactic longitude to Sgr, and towards the south- 
ern extension of Sgr, its motion towards the Sgr, and its orbital 
parameters i ndicate that it may ha ve been captured by Sgr some 
1.7 Gyr ago (iDinescu et al.Lll997h . 

From the chemical abundance point of view, a few 

S Fe/H1 e stimates place it in the m oderately metal-poor regime 
Brown et al, 19 97; Cohen. |2004 see also Table [3]l. In particu- 
lar, CohenI ( 20041) noticed not only the lack of a-enhancement in 
Pal 12, but also remarked on the abunda nce of heavy elements, 
which is more similar to Sgr dwarf s tars dBonifacio et all 120001 : 
ISmecker-Hane & M cWilliaml l2002l) . The four stars she studied 
showed identical abu ndance pattern, with no evidence of star-to- 
star variations. Also iKavser et al.l (l2008h found no sign of CH 
and CN anti-correlation among evolved stars. 

With our data, a bimodal anti-correlation of CH and CN band 
strengths is apparent. The group of CH-strong and CN-weak 
stars appears truncated and reduced in size, probably owing to 
our selection (Sect. [23T l. which excluded all stars with S/N<8 in 
the CN spectral region. A higher S/N sample could shed some 
light on this peculiar effect, which makes the CN-strong to CN- 
weak stars ratio appear quite higher than in other clusters (see 
also Sect.|6|l. 



6. Trends with cluster parameters 

When trying to connect the GGC properties with the presence 
and extension of the anti-correlations, two parameters are usu- 
ally defined: (i) the "low -resolution community", studying the 
strength of molecular indices, u sually builds the ratio of CN- 
strong to CN-weak stars Tcn (e.g.lNorrisl[l987l;ISmith & Mateol 



[199 



SmithLl2002HKayser et alil2008h ; (ii) the "high-resolution 
community" (namelv ICarrettal l2006t ICarretta et al]. l2007cl) 



measures the extension of the Na-O and Mg-Al anti-correlations 
with interquartile ranges. The two parameters measure two dif- 
ferent physical quantities, and both give clues to the understand- 
ing of the anti-correlation phenomenon. In this paper we will 
concentrate on the first approach, deferring the second one to a 
following paper, in which we will attempt to derive [C/Fe] and 
[N/Fe] ratios through spectral synthesis of the best S/N clusters. 



Cluster 






ICN 


Stcn 


ACN 


SACN 


NGC104 


33 


31 


1.06 


0.27 


0.209 


0.051 


NGC288 


60 


69 


0.87 


0.15 


0.088 
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To compare with Tcn, we selected several cluster intrinsic 
(Table O and orbital parameters (Table HJi from the literature. 
Because we were looking for trends and were not particularly 
interested in zeropoints and absolute quantities, we preferred pa- 
pers reporting homogeneous data, trying to avoid mixing of dif- 
ferent sources whenever possible. 

6.1. Ratio of CN-strong to CN-weak stars 

The ratio between CN-strong and CN-weak stars tells us the rela- 
tive importance of stellar groups with different chemical compo- 
sition. In the framework of the self-enrichment scenario, for ex- 
ample, it becomes a fundamental constraint to model the chem- 
ical evolution of a GGC. A correlation between Icn a nd clus- 
ter el li pticity, defined a s e-l-(b/a), was init ially found jNorrisl 
[19871 lsmith & Mateol [l990l Isiniflil l200l. but later studies 
('Kayser et al.', '2008*) did not confirm it. ISmith & Mateol (Il990h 
and Kayser et al. (2008) found that clusters with a higher frac- 
tion of CN-strong stars are also more luminous and therefore 
more massive. T his finding is supported by the trend found by 
ISmith & Mateo (1199 0) with the central velocity dispersion ctq. 
While iKavser et al.l (l2008h . in their compilation of literature es- 
ti mates of Tcn, found al l po ssible values o f Icn be tween and 
2, ICarretta etd] ( l2009ah and lCarretta et al.l (l20()9bh . using more 
homogeneus data, found that the second - or I as intermediate - 
generation stars (i.e., the CN-strong stars) are between 50% and 
70% of the total. 

We computed our Tcn ratio by selecting stars with 
5S3839(CN)>0 for our CN-strong group and stars with 
5S3839(CN)<0 for our CN-weak group. The final errors on Ycn 
were propagated by assuming Poissonian errors ( y/n) in the CN- 
strong and CN-weak count£| and are reported in Table |5] We 
found that even if the CH4300 region has a much higher S/N in 
our spectra, the spread of a corresponding Ich ratio is too large 
to identify any correlation. We interpret this as a consequence of 
the fact that CN, being a double-metal band, is more sensitive to 
abu ndance variati ons and therefore Tcn is the best indicator (see 
also lNorrisL[T987h . 



d = s/lx^ -I- y^), where x is the difference between the average 
^S3839(CN) of the (5S3839(CN)>0 and (5S3839(CN)<0 groups, while 
y is the difference between the average 5CH4300 of the <5S3839(CN)>0 
and (5S3839(CN)<0 groups. The two centroids in the CN-CH plane are 
represented by the two lage dots in each panel of Fig.jTj ACN is reported 
along with its uncertainty in the last two columns of TablejS] 



' As noted in Sect. 13.21 the very small uncertainty in the placement 
of the sparation line between CN-strong and CN-weak stars can have 
an impact on Icn in the sense of moving a few stars from one group to 
the other. This uncertainty ranges from the order of 1% to a maximum 
of 10%. In any case, it is always small compared to the uncertainties 
quoted in Table [5] 
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Fig. 10. Run of the ratio TcN-'ncN„„,„J'^CN,,„,t for our sample clusters. In each panel, NGC 362 ans Pal 12 (grey dots) are excluded 
from the relations. Dotted lines mark linear fits to the remaining clusters (black dots) and the correlation coefficients are reported on 
top of each panel: r^ stands for the Spearman and rp for the Pearson coefficient. For the meaning of X labels, see Tables[3]and|4l 
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The average ratio we find between the two populations, 
<rcA'>= 0.82+0.29, implies that the CN-strong stars - the sec- 
ond generation - tend to constitute the minority of the total po p- 
ulation, in contrast with the results by ICarretta et al.l (l2009allbh . 
In both studies there are ample uncertainties, with our esti- 
mated fraction of second generation stars varying from roughly 
25% to 50% of the total, and the same fraction estimated by 
ICarretta etall ( 2009a b) varying from about 50% to 70% of the 
total. While the discrepancy is only marginal, and indeed the 
data appear consistent within the uncertainties, it is nevertheless 
puzzling that these two indicators - the CH-CN and the Na-O 
anti-correlations - give somewhat different answers (see also the 
discussion in Sect. 14. U . 

6.2. Comparison with cluster parameters 

Our results are reported in Fig.[TOl where we also report on the 
Spearman (r^) and Pearson (r/>) correlation coefficient^ As a 
first note, we observe that our Tcn values all range between 0.5 
and 1.0 (<rcAf>=0.82±0.29), with the only exceptions of Pal 12 
and NGC 362, which have much higher uncertainties than the 
rest of the clusters. Indeed, our strict selection of spectra with 
S/N>8 could have biased our results for Pal 12 against CN-weak 
stars, while for NGC 362 the S/N ratio and the low number of 
stars prevent us from measuring rcAf properly (see also the error- 
bars in Fig.|6]l. Therefore, we have decided to leave out NGC 362 
and Pal 12 from the present attempt to look for linear correlations 
of Tcn and cluster parameters. 

We do not find any highly significant correlation of clus- 
ter parameters with Icn- The only correlations that have sig- 
nificance slightly above 1 cr, are the ones of Icn with cluster 
age (83%, or ^1.2cr), concentration c (88%, or ^1.3cr), total 
luminosity My (82%, or ^1.2 o"), total present-day mass in units 
of log(M/Mo) (81%, or ^1.2 cr), orbital inclination (95%, or 
-2(t), and orbital period P (78%, or ^1.1 cr). However weak, 
these relations broadly agree with past findings, except for the 
correlation between Tcn and ellipticity. In the framework of the 
self-enrichment scenario, this might mean that the relative size 
of the two stellar generations (or the efficiency of the second star 
formation episode) is only marginally influenced by the cluster 
total mass and orbital parameters. The correlation with cluster 
age is detected here for the first time and if confirmed, it could 
mean that older clusters are more efficient in when converting 
gas into the second generation of stars (but see Sect. |7]i. 

7. Discussion and conclusions 

We can summarize the results of this paper as follows: 

- for the metal-rich clusters (Pal 12, 47 Tuc, NGC 5927, 
NGC 6352) where the S/N of our spectra was good enough, 
we could detect clear anti-correlations between the CN and 
CH band strengths, which were always clearly bimodal; 

- for the metal-poor clusters we could detect no clear anti- 
correlation or bimodality, even in those cases where we knew 
from past studies that variations should be present; we as- 
cribe that to the fact that we were unfortunately not able to 
gather sufficient S/N (although higher than for the metal-rich 
clusters) to study CH and CN bands in these metal-poor MS 
stars; 

* To compute the Spearman and Pearson coefficients a nd their as- 
sociate d probabilities, we used the algorithms described bv lPress et alj 
(Il992h in their "Numerical recipes". 



- the only exception was M 15, the most metal-poor of the 
entire sample; which shows an abnormally large difference 
in band strengths for its metallicity; 

- we computed Icn, the ratio between the CN-strong and 
CN-weak stars, which turned out to be on average 
<TcN >=0.82+0.29; 

- no strongly significant correlation was found between Icn 
and 15 clusters parameters deduced from literature compila- 
tions; 

- weak correlations (slightly above 1 cr) were found with clus- 
ter total present-day mass and luminosity; 

- similarly weak correlations were found with orbital parame- 
ters, with cluster concentration and with cluster age. 

Now, as reviewed in Sect. [1] the so-called intrinsic scenar- 
ios involving various kinds of extra-mixing are not sufficient to 
explain the observed chemical anomaUes in GGC, although o f 
course mixing does occur in GGC stars (iGratton et al.L l2000l) . 
This is especially true if one considers, as we did, only un- 
evolved stars that still have not undergone the so-called first 
dredge-up episode. We will therefore focus on the extrinsic sce- 
narios, and in particular on the self- enrichment scenario. 

Our first result is that in all the studied clusters, those which 
have sufficiently good data to show CH and CN anti-correlations 
among MS stars are clearly bimodal in both the CN and CH 
band strengths. This gives a very strong support to the above 
scenario, especially when it predicts two (or more) subsequent 
stellar generations, each polluting the gas from which the fol- 
lowing one forms. In this respect we can say that within our 
sample, we always find two distinct groups of stars and not more. 
So, while exceptions are always possible (let us not forget about, 
e.g., (xi Cen and possibly NGC 2808), the general tendency seems 
to have at most two major stellar generations in each GGC0. 

Our second result is that M 15 appears to possess exception- 
ally large variations of the CH and CN band strengths, for its 
metallicity. If this, as expected, will lead to an exceptionally ex- 
tended [C/Fe] vs [N/Fe] anti-correlation, it would certainly be 
worthwhile to follow-up on M 15 with large samples of high- 
resolution spectra, to better study the presence of multiple pop- 
ulations also in the Na-O plane. Given its peculiar properties 
(Sect. 15.111 ). M 15 could even belong to the group of massive 
clusters that have undergone more self-enrichment than usual, 
either because of their special properties or because they could 
be the remnants of accreted dwarf galaxies (such as o) Cen). 
Therefore, also further searches for variations in heavier species 
such as Ca, Fe, or Ba should be extremely interesting. 

Our third result is that the ratio of CN-strong to CN-weak 
stars is on avera ge <rcN >-0.82±0.29. T his challenges to some 
extent results bv lCarretta et alj ( l2009albl) . who find that the first 
or primordial (P) generation accounts for approximately 1/3 of 
the total cluster population, leaving space for a 50-70% fraction 
for the second or intermediate (1) population. On the one hand 
the discrepancy is - at best - marginal, as discussed in Sect. |6] 
and given the ample uncertainties involved; on the other hand, 
the discrepancy could be caused by the Na-O and the CH-CN 
anti-correlations sampling two sligthly different phenomena (as 
discussed in Sect. 14. 11 1. Certaintly more studies are required to 
shed more light on this issue, because these kind of discrepancies 



' Of course, we cannot eclude with our study the existence of small 
populations totalling less than ^10-20% of the respective cluster pop- 
ulations. Therefore, t he existance of an extreme (E in the notation by 
ICarretta et all l2009al lbl) population is by no means excluded by our 
data. 
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- if confirmed - could turn out into further clues on the formation 
and evolution of GGC. 

Our fourth result is that no strong correlation was found be- 
tween rcN ond 15 different cluster parameters compiled from the 
literature. While we must stress that the present study is limited 
to 12 clusters with spectra of modest S/N, the absence of strong 
corelations has one major implication in the self- enrichment sce- 
nario: it would imply that the population ratio of the two (or 
more) stellar generations is only mildly (if at all) influenced by 
the cluster properties. In particular, clusters with larger present- 
day total masses (and therefore luminosities) could be able to 
form only sligthly higher fractions of second generation stars 
with respect to less massive clusters. Similarly, the relations with 
cluster orbital parameters could hint that environment has some 
(minor) role in determining the size of the second stellar gener- 
ation. 

In summary, the use of low-resolution spectra to study CH 
and CN bands in GGC stars provides a useful complement to 
the high-resolution studies, needing less observing time and data 
analysis effort. In particular, the use of strictly unevolved stars, 
i.e., stars that have not undergone first dredge-up mixing, and 
of two elements such as C and N that (unlikely Na and O) 
are strictly produced in the same nuclear reaction chains and at 
the same temperature, can give additional constraints to studies 
based on red giants only. 
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